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Abstract 



o 
o 
o 

■ The simplest unified extension of the Minimal Supersymmetric Standard Model 

with bi-linear R-Parity violation naturally predicts a hierarchical neutrino mass 
spectrum, in which one neutrino acquires mass by mixing with neutralinos, while 
the other two get mass radiatively. We have performed a full one-loop calculation of 
the neutralino-neutrino mass matrix in the bidinear ^ P MSSM, taking special care to 
^ ' achieve a manifestly gauge invariant calculation. Moreover we have performed the 

renormalization of the heaviest neutrino, needed in order to get meaningful results. 
The atmospheric mass scale and maximal mixing angle arise from treedevel physics, 
while solar neutrino scale and oscillations follow from calculable one-loop correc- 
tions. If universal supergravity assumptions are made on the soft-supersymmetry 
breaking terms then the atmospheric scale is calculable as a function of a single 
Ijip violating parameter by the renormalization group evolution due to the non-zero 
bottom quark Yukawa coupling. The solar neutrino problem must be accounted for 
by the small mixing angle (SMA) MSW solution. If these assumptions are relaxed 
then one can implement large mixing angle solutions, either MSW or just-so. The 
theory predicts the lightest supersymmetic particle (LSP) decay to be observable at 
high-energy colliders, despite the smallness of neutrino masses indicated by exper- 
iment. This provides an independent way to test this solution of the atmospheric 
and solar neutrino anomalies. 



1 Introduction 



The high statistics data by the SuperKamiokande collaboration [Q has confirmed the 
deficit of atmospheric muon neutrinos, especially at small zenith angles, opening a new 
era in neutrino physics. On the other hand the persistent disagreement between solar neu- 
trino data and theoretical expectations has been a long-standing problem in physics @. 
Altogether these constitute the only solid evidence we now have in favour of physics 
beyond the present standard model, providing a strong hint for neutrino conversion. Al- 
though massless neutrino conversions || can be sizeable in matter, and may even provide 
alternative solutions of the neutrino anomalies [Q, it is fair to say that the simplest 
interpretation of the present data is in terms of massive neutrino oscillations. Taking 
for granted such an interpretation, the present data do provide an important clue on the 
pattern of neutrino masses and mixing. The atmospheric data indicate to v T flavour 
oscillations with maximal mixing ||, while the solar data can be accounted for in terms 
of either small (SMA) and large (LMA) mixing MSW solutions |J, as well as through 
vacuum or just-so solutions [|7|]. A large mixing among v T and v e is excluded both by 
the atmospheric data and by reactor data on neutrino oscillations ||. There has indeed 
been an avalanche |§ of papers trying to address this issue in the framework of unified 
models adopting ad hoc texture structures for the Yukawa couplings. 

Here we propose an alternative approach to describe the structure of lepton mixing 



which accounts for the atmospheric and solar neutrino anomalies [|10| based on the simplest 
extension of minimal supergravity with bi-linear R-Parity violation JTI ]. The particles 
underlying the mechanism of neutrino mass generation are the neutral supersymmetric 
partners of the Standard Model gauge and Higgs bosons which have mass at the weak-scale 
and are thus accessible to accelerators. 

Our model breaks lepton number and therefore necessarily generates non-zero Ma- 
jorana neutrino masses fl2fl . At tree- level only one of the neutrinos picks up a mass by 
mixing with neutralinos JT3 |, leaving the other two neutrinos massless ]H| . While this 



can explain the atmospheric neutrino problem, to reconcile it with the solar neutrino data 
requires going beyond the tree-level approximation. This is the purpose of the present pa- 
per. Here we improve the work of ref. |IJ| by performing a full one-loop calculation of the 
neutrino mass matrix and also update the discussion in the light of the recent global fits 
of solar and atmospheric neutrino data. This can also be used to improve the discussion 
given in |]TB| where the tree approximation was assumed. For simplified analyses including 



only the atmospheric neutrino problem in the tree-level approximation see ref. |T7j and a 



number of papers in ref. [[18 



We have performed a full one-loop calculation of the neutralino-neutrino mass matrix 
in the bi-linear ]fl p MSSM, showing that, in order to explain the solar and atmospheric 
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neutrino data, it is necessary and sufficient to work at the one-loop level, provided one 
performs the renormalization of the heaviest neutrino. In contrast to all existing papers 
|T5| , jlBf , we have taken special care to verify the gauge invariance of the calculation, thus 
refining the approximate approaches so far used in the literature. We find that if the 
soft-supersymmetry breaking terms are universal at the unification scale then only the 
small mixing angle (SMA) MSW solution to the solar neutrino problem exists. On the 
other hand if these assumptions are relaxed then one can implement large mixing angle 
solutions, either MSW or just-so. 

Bilinear R-parity breaking supersymmetry has been extensively discussed in the 



literature ||10|| . It is motivated on the one hand by the fact that it provides an effec- 
tive truncation of models where R-parity breaks spontaneously by singlet sneutrino vevs 
around the weak scale |19| |. Moreover, they allow for the radiative breaking of R-parity, 
opening also new ways to unify Gauge and Yukawa couplings p0| and with a potentially 
slightly lower prediction for a s ||21|| . For recent papers on phenomenological implications 
of these models see ref. |22| , p3| , If present at the fundamental level tri-linear breaking 
of R-parity will always imply bi-linear breaking at some level, as a result of the renor- 
malization group evolution. In contrast, bi-linear breaking may exist in the absence of 
tri-linear, as would be the case if it arises spontaneously. 

This paper is organized as follows. In sections 0, [3] and [| we describe the model, 
the minimization of the scalar potential and the radiative breaking of the electroweak 
symmetry. In section [5] the tree level masses and mixings are described, while the contri- 
butions to the one loop mass matrix and the gauge invariance issue are studied in section 
H Finally the neutrino masses and mixings are discussed in section [7] where we show our 
results for solar and atmospheric oscillation parameters. The more technical questions 
regarding the mass matrices, couplings and one loop results as well as further details of 
gauge invariance are given in the appendices. We also briefly discuss how, despite the 
smallness of neutrino masses indicated by experiment, the theory can lead to observable 
I/ip phenomena at high-energy accelerators. 



2 The Superpotential and the Soft Breaking Terms 



Using the conventions of refs. |23|, |25| we introduce the model by specifying the superpo- 
tential, which includes BRpV |IU| in three generations. It is given by 



W = e ab h^QIUfil + h%Q h j5jE a d + h%t\R 3 H a d - ^H a d H b u + e^H b u 



where the couplings hu, hp and He are 3x3 Yukawa matrices and /i and q are parameters 
with units of mass. The bilinear term in eq. (|l|) violates lepton number in addition to 
R-Parity. 
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Supersymmetry breaking is parameterized with a set of soft supersymmetry breaking 
terms. In the MSSM these are given by 

l mssm = M^QfQ* + Mfufi] + MfDiD* + Mflfl) + Mf^R) 

+m 2 H H d l *H a d + m 2 u H a *H a u - \\M S \ S \ S + \M\\ + \M'\'\' + h.c] (2) 



+£ab 



A^UjH b u + A%Q\D 3 H a d + A%LiR,H$ - B^H a d H b u 



In addition to the MSSM soft SUSY breaking terms in C^f the BRpV model contains 
the following extra term 

VZT = -BieieJ.tK , (3) 

where the Bi have units of mass. In what follows, we neglect intergenerational mixing in 
the soft terms in eq. (H). 

The electroweak symmetry is broken when the two Higgs doublets Hd and H u , and 
the neutral component of the slepton doublets L\ acquire vacuum expectation values. We 
introduce the notation: 

-(3)' *•-(%)• Mf)< 

where we shift the neutral fields with non-zero vevs as 

H^^K + va + iy^, H u = -L[a u + v u +itp u }, L° = ±={p* + v, . (5) 



Note that the W boson acquires a mass = \g 2 v 2 , where v 2 = v 2 d + v 2 + v \ + v f + v 2 ~- 

\2 



(246 GeV) . We introduce the following notation in spherical coordinates for the vacuum 
expectation values: 

v d = v sin 6*i sin 2 sin 9 3 cos ft 
v u = v sin 9\ sin 9 2 sin 63 sin ft 

v 3 = v sin #1 sin 9 2 cos #3 (6) 
v 2 = v sin 6*i cos 6 2 
Vi = v cos 61 

which preserves the MSSM definition tan/3 = v u /vd- In the MSSM limit, where = 
Vi = 0, the angles 9{ are equal to tt/2. In addition to the above MSSM parameters, our 
model contains nine new parameters, e*, Vi and Bi. The three vevs are determined by the 
one-loop tadpole equations, and we will assume universality of the 5-terms, B = Bi at 
the unification scale. Therefore, the only new and free parameters can be chosen as the 
e, ; . 
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3 The Scalar Potential 



The electroweak symmetry is broken when the Higgs and lepton fields acquire non-zero 
vevs. These are calculated via the minimization of the effective potential or, in the 
diagramatic method, via the tadpole equations. The full scalar potential at tree level is 

2 



v° = V 

v total / j 



dW 



dzi 



+ Vn + Cf M + VfoT (7) 



where Zi is any one of the scalar fields in the superpotential in eq. ([!]), Vd are the D-terms, 
and V^ff is given in eq. (|3[). 

The tree level scalar potential contains the following linear terms 

^linear = td^d + tu^u + h v \ + h v 2 + > ( 8 ) 

where the different t° are the tadpoles at tree level. They are given by 



*d = ( m L + ^) V d + V d D - h(Bv u + Vi€i 

£° = -B/2v d + (m 2 Hu + /j 2 ^v u -v u D + v i B i e i + v u e 2 

tl = v x D + e x ( - [iv A + v u Bx + v^) + \ (viM 2 Lil + M| H ^) (9) 



v, 



t° 2 = v 2 D + e 2 (- nv d + v u B 2 + V&) + \ (viM 2 Li2 + M- 
4 = v 3 D + e 3 (- fiv d + v u B 3 + v^) + \ (viM 2 Li3 + M 2 L3i 

where we have defined D = \(g 2 + g' 2 ){v\ + v\ + v\ + v\ — v 2 ) and e 2 = e\ + e\ + e\. A 
repeated index i in eq. (|9|) implies summation over i = 1,2, 3. The five tree level tadpoles 
t° a are equal to zero at the minimum of the tree level potential, and from there one can 
determine the five tree level vacuum expectation values. 

It is well known that in order to find reliable results for the electroweak symmetry 
breaking it is necessary to include the one-loop radiative corrections. The full scalar 
potential at one loop level, called effective potential, is 

Vtotal = V t ° otal + V RC (10) 

where Vrc include the quantum corrections. In this paper we use the diagramatic method, 
which incorporates the radiative corrections through the one-loop corrected tadpole equa- 
tions. The one loop tadpoles are 

t a = t°a ~ + UQ) =t° a + 7™{Q) (11) 

where a = d,u, 1,2,3 and T® R (Q) = —6t^ IS + T a (Q) are the finite one loop tadpoles. 
At the minimum of the potential we have t a = 0, and the vevs calculated from these 
equations are the renormalized vevs. 
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Neglecting intergenerational mixing in the soft masses, the five tadpole equations 
can be conveniently written in matrix form as 



j.0 + j.0 j.0 j.0 
L ui b di L li h 2i L 3 



where the matrix M 2 ad is given by 

, 2 



M t 2 ad 



M tad [Vu,V d ,V 1 ,V 2 ,V 3 } J 



;i2) 



d +^ + D 


-Bfl 


-/iei 


~l^e 2 




-Bfi 




-D B lGl 


B 2 € 2 


B 3 e 3 




Sid 


Ml+ej+D 


eie 2 


eie 3 




B 2 € 2 


eie 2 


Ml 2 +el+D 


C2C3 


-/^ 3 


B 3 e 3 


ei£3 


e 2 e 3 


M 2 L3 +el + D_ 



(13) 



and depends on the vevs only through the D term defined above. 

In order to have approximate solutions for the tree level vevs, consider the following 
rotation among the H d and lepton superfields: 

1,2 RM t 2 ad R 1 



M' 



tad 



(14) 



where the rotation R can be split as 
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1. 


where the three 


ang 


les 


are defined as 




ei 
A*" 




























Cl = — , 

n' 




Si = 






= V 


'/*' 


+ Cl, 


















c 2 = — , 

n" 




s 2 = 


e 2 
V? 




A*" 


= \ 


/V 2 + el 









(15) 



(16) 



c 3 



A* e 3 



/i" 2 + e|. 



It is clear that this rotation R leaves the D term invariant. The rotated vevs are given 
by 



l V L V 'd, V 'li V 2, v'sl = R K> ^ v i, v 2, v 3 f , 



(17) 



and under the assumption that v' 1 ,v' 2 , v' 3 <C v, these three small vevs have the approximate 
solution 



Vn ~ 



m 2 — Mr , 



/x'/x'" 



v d + 



A* e 2 



A^ 3 
Mg+D 



m 



Ml . 
v d + 



n"n"' 



v 



/// 



///2 



d 



B x -B , 

B 2 -B' , 

B 3 - B" 
v 

fJL'" 



(18) 
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where we have defined the following rotated soft terms: 

n _ rn^y + MlA ,,2 _ + MIA m m _ + M IA 

lib TT "~ . I I L TT ~ . Ill 



"H d - n > " b H d - „ 2 ' "°H d - m > 

By? + B x t\ g> /2 + B 2 e\ By" 2 + B 3 e 2 
B = ^ . B = , B = ^ . (19) 

The approximation v'^v'^v^ <C v is justified in SUGRA models with universality of soft 
masses at the weak scale, as shown in the next section. 



4 Radiative Breaking of the Electroweak Symmetry 



It was demonstrated in ref. JT(J that BRpV can be succesfully embedded into SUGRA 
with universal boundary conditions at the unification scale, and with a radiatively bro- 
ken electroweak symmetry. At Q — Mjj we assume the standard minimal supergravity 
unification assumptions, 

A t = A b = A T = A , 
B = B i = A - 1 , 

m\ d = m 2 Hu = Ml = Ml = Mq. = My. = M 2 Di = m\ , (20) 
M A M, .I/, Mi/a • 

We run the RGE's from the unification scale My ~ 2 x 10 16 GeV down to the weak scale, 
giving random values to the fundamental parameters at the unification scale: 

1(T 2 < h\ v I^K < 1 

< h 2 bu /iir < 1 

-3 < a = A/m < 3 (21) 

< nfr/m 2 < 10 

< M 1/2 /m < 5 

The Yukawa couplings are determined by requiring that three eigenvalues of the chargino / charged- 
lepton mass matrix corrrespond to the experimentally measured tau, muon, and electron 
masses f\ 

As in the MSSM, the electroweak symmetry is broken because the large value of 
the top quark mass drives the Higgs mass parameter m 2 Hu to negative values at the weak 

*For the case of large tree-level neutrino mass one must note that the lepton Yukawa couplings are no 
longer related to the lepton masses via the simple relations valid in the Standard Model. Since charginos 
mix with charged leptons, the Yukawa couplings depend also on the parameters of the chargino sector. For 
the case of interest here (light v r mass fixed by the atmospheric scale) this correction is less important. 
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scale via its RGE [26]. In the rotated basis, the parameter fi" 12 is determined at one loop 

by 

1 r ~ i ( m 'u + T5 R )-(m 2 + TS R ) t' 2 
V ~ ~2 l z ~ Azz ( m z)\ H 1 / 2 _ l 

where t'g = v' u /v' d is defined in the rotated basis and is analogous to tan/3 in eq. (||) 
defined in the original basis. The finite DR Z-boson self energy is Azz(m|), and the 
one-loop tadpoles T® R and T^ R are obtained by applying to the original tadpoles in 

d u 

eq. fllT]) the rotation R defined in eq. ((0^). The radiative breaking of the electroweak 
symmetry is valid in the BRpV model in the usual way: the large value of the top quark 
Yukawa coupling drives the parameter m 2 Hu to negative values, breaking the symmetry of 
the scalar potential. 

As we will see a radiative mechanism is also responsible for the smallness of the 
neutrino masses in models with universality of soft mass parameters at the unification 
scale. The relevant parameters are the bilinear mass parameters B and Bi, the Higgs 
mass parameter rn 2 Hd , and the slepton mass parameters M 2 .. 

The RGE's for the B parameters are 

^ = JL (3h 2 A t + 3h 2 b A b + h 2 T A T + 3g\M 2 + f^Mx) 

it = i + k " Ar + 3 & M * + ^ ?Ml ) (23) 

dB 2 



^ = -L (3h 2 A t + 3g\M 2 + lg 2 Ml 



dt dt 8tt 2 

where we do not write the effect of Yukawa couplings of the first two generations. Similarly, 
the RGE for the down-type Higgs mass is 

^ = JL (3h 2 X b + h 2 T X T - 3g 2 2 M 2 - lg 2 M 2 ) , (24) 
and the RGE's for the slepton mass parameters are 

2 



- % {3g\Ml + IglM 2 ) , (25) 



dt dt 8tt 2 

where X 6 = m 2 Hd + Mj 3 + M 2 Di + A 2 and X r = m 2 Hd + M| 3 + M| 3 + A 2 . 

With the aid of these RGE's we can find an approximate expression for the slepton 
vev's in the rotated basis v' t , given in eq. flTS|). The relevant soft term differences, defined 
as ABi = Bi — B and Am- = M| . — m 2 ^ , are approximated by 



Afi 3 = -^(3^)ln^ 



A£? 2 = AB 1 = -L (3/i 6 2 A 6 + ^ A r ) In (26) 
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for the B terms, and by 

Mi 



8vr 2 V i m weak 



Ami = Am 2 = -L (3h 2 X b + h 2 T X T ) In (27) 
8tt 2 v / m weak 

for the mass squared terms. This way, if we assume that <C /i we can neglect the 
rotations in eq. (|l8|) and we find 



Ml +D 



Am- — tpfxABA (28) 



which give us an approximate expression for the sneutrino vev's v[ in the basis where the 
Cj terms are absent from the superpotential. In a model with unified universal boundary 
conditions on the soft SUSY breaking terms (SUGRA case, for short) the v\ are calculable 
in terms of the renormalization group evolution due to the non-zero bottom quark Yukawa 
coupling. We should stress here that for our subsequent numerical calculation we solve 
the tadpole equations exactly. 

The symmetry of the neutralino/neutrino mass matrix implies that only one neutrino 



acquires a tree level mass, and the other two remain massless jnj (see next section). The 
massive neutrino will have the largest component along r, /i or e if the largest vev is v' 3 , 
v' 2 , or v[ respectively. On the other hand, the most obvious difference between the third 
generation sneutrino vev and the first two generations is in the extra contribution from 
h T to ABi in eq. ( ^6l) and to Am- in eq. ( ^7|) for the first two generations. Due to the 
tau lepton contribution, AB X and AB 2 are larger than AB 3 , and similarly for the Am-, 
specially if tan/3 ^> 1. However, we have checked that it is possible without fine-tuning 
the parameters in an unnatural way to arrange for the heaviest of the neutrinos to be an 
equal mixture of and u T as needed in order to obtain an explanation of the atmospheric 
neutrino anomaly. That this is possible can be understood by noticing that there can be 
a cancellation between the AB and Am 2 terms in eq. (|28|) for v[ and v' 2 . 



5 Tree Level Neutrino Masses and Mixings 

Here we discuss the tree level structure of neutrino masses and mixings. For a complete 
discussion of the fermion mass matrices in this model see Appendix [A|. f\ In the basis 

4n our notation the four component Majorana neutral fermions are obtained from the two component 
via the relation 
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or 



—i\ 3 , Hj, Vei v t) the neutral fermion mass matrix M N is given by 

M N = ' (29) 



m 







where 



Mi 








M 2 

-\gv u 



\gvd 



\g'v u 
-\gv u 



(30) 



is the standard MSSM neutralino mass matrix and 

-\g' v i \g v i o e i 

- \g'v 2 \gv 2 e 2 

- iafva \gv-i e 3 



m 



(31) 



characterizes the breaking of R-parity. The mass matrix M\- is diagonalized by (see 
Appendix |A|) 

N*M N N~ l = diag(m x o, m „ ) (32) 



where (i = 1, • • • , 4) for the neutralinos, and (j = 1 



3) for the neutrinos. 



We are interested in the case where the neutrino mass which is determined at the 
tree level is small, since it will be determined in order to account for the atmospheric 
neutrino anomaly. The above form for is especially convenient in this case in order 
to provide an approximate analytical discussion valid in the limit of small Tjl p violation 
parameters. Indeed in this case we perform a perturbative diagonalization of the neutral 
mass matrix, using the method of p7f, by defining p4j 



i = m ■ Mj 



If the elements of this matrix satisfy 



« 1 



(33) 



(34) 



then one can use it as expansion parameter in order to find an approximate solution for 
the mixing matrix N '. Explicitly we have 



6i 

£i3 



g'M 2f i 

2det(M x n 

gM lfJ , 
2det(M x o 



■A; 



-A. 



-- + 



(g 2 M 1 +g' z M 2 )v u 
Adet(M x o) 



A; 
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(g Mi + 9 M 2 )vd 
4det(M x o) 



A; 



where 



Aj = fIVi + V d €i oc vt 



(35) 



(36) 



are the alignment parameters. From eq. (^) and eq. (^) one can see that £ = in the 
MSSM limit where = 0, Vi = 0. In leading order in f the mixing matrix M is given by, 



Af*-- 



N* W 1 - if tf ft 

o K T ) { -f l - kft 



(37) 



The second matrix above block-diagonalizes the mass matrix M n approximately to the 
form diag(.M x o, m e //), where 



-m ■ Ai« Tn 

X, 



M l9 2 + M 2 g 
4det(M x o) 



( A 2 e A e A, A e A T \ 
Af. A„A T 



A e A M 



(3? 



\A e A T A M A r A^ J 



The sub-matrices A/" and V v diagonalize -M x o and m e // 

N*M x oN^ = diag(m x o), (39) 
V?m eff V v = diag(0, 0, m„), (40) 

where 

M ig 2 + M 2 g' 2 ltl2 
m v =Tr(m eff ) = \M • (41) 

Clearly, one neutrino acquires mass due to the projective nature of the effective neutrino 
mass matrix m e ff, a feature often encountered in ]fl p models [14]. As a result one can 
rotate away one of the three angles |I2] in the matrix V„, leading to [^ 

/ 1 



v v = 







\ / cos 6»i3 -sin 6»i3 \ 
x 1 



cos 6*23 — sin 9 2 % . 
\ sin #23 cos #23 / \ sin 13 cos 



#13 / 



(42) 



where the mixing angles can be expressed in terms of the alignment vector A as follows: 

A, 



tan# 



13 



tan 9 



23 



At 



(43) 
(44) 
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6 One Loop Neutrino Mass Matrix 



One-loop radiative corrections to the neutralino/neutrino mass matrix in the BRpV model 
were calculated first in [IS], working in the t'Hooft-Feynman gauge (£ = 1). Our analysis 
improves the previous work in that we check explicitly the gauge invariance using the 
gauge. We use dimensional reduction to regularize the divergences |29| and include all 
possible MSSM particles with consistently determined mass spectra and couplings in the 
relevant loops. 



6.1 Two— Point Function Renormalization 



We denote the sum of all one-loop graphs contributing to the 2-point function as 




The most general expression for the one-loop contribution to the unrenormalized neu- 
tralino/neutrino two-point function is 



p l tiUp 2 ) + P R n?Ap 2 



(45) 



where the indices i and j run from 1 to 7, Pr = |(1 + 75) and Pl = |(1 — 75) are the 
right and left projection operators, and p is the external four momenta. The functions S 
and I! are unrenormalized self energies and depend on the external momenta squared, p 2 . 
The neutral fermions Ff are a mixture of weak eigenstate neutralinos and neutrinos and 
given by 

^° = A/^° (46) 

where M is the 7x7 matrix that diagonalizes the neutralino/neutrino mass matrix ac- 
cording to eq. (|32|) . 

The inverse propagator at one loop is obtained by adding to the tree level propagator, 



this self energy previously renormalized with the dimensional reduction DR scheme and 



denoted as E and II. In the DR scheme, the counterterms cancel only the divergent pieces 
of the self energies. In this way, they become finite and dependent on the arbitrary scale 
Q. The tree level masses are promoted to running masses in order to cancel the explicit 
scale dependence of the self energies. Thus, the inverse propagator of the neutral fermion 
is 



F? 



?Pf(p) = ~ ™fM + S^(p, Q) 



(47) 
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The physical pole mass is given by the zero of the inverse propagator, in the limit where 



m Fi , and may be found using 



Z F t U (P) 



u(p) = u{p) prf - m Fi (Q) + Z FF (p, Q) u(p) 



(48 



where u and u are two on-shell spinors, m F . and m Fi (Q) are the neutral fermion pole and 
running masses respectively, and ~S FF (p, Q) is the renormalized two-point function in the 
DR scheme. The quantity Z F corresponds to the finite ratio of the infinite wave function 
renormalization constants in the DR scheme and the on-shell scheme, and it accounts 



for the fact that the residue of the DR propagator at the pole is not one ||30|| . The 
renormalized function ~E FF (p, Q) is calculated by subtracting the pole terms proportional 
to the regulator of dimensional reduction 

2 



4-d 



7s + ln47r 



(49) 



where j F is the Euler's constant and d is the number of space-time dimensions. In practice 
we have 



T$ F (p,Q)= ?? FF (p) 



A=0 



(50) 



Since u{p)^u{p) 
eq. (ESI) we find 



= 0, the terms proportional to 75 in "S FF do not contribute. From 

Am Ft = m Fi - m Fi (Q) = f^(mi) - m Fi S^(m|.) , (51) 



where 



n 



(52) 



and the tilde implies renormalized self energies. A given set of input parameters in the 
neutralino/neutrino mass matrix defines the set of tree level running masses m F .(Q), 
among them two massless and degenerate neutrinos. The one-loop renormalized masses 
m Fi are then found through eq. (|5T|), and the masslessness and degeneracy of the two 
lightest neutrinos is lifted. 

The tree level masslessness of the lightest neutrinos implies an indetermination of 
the corresponding eigenvectors. In order to find the correct neutrino mixing angles we 
diagonalize the one-loop corrected neutralino/neutrino mass matrix. We define 



Mi 



pole 



Mg R (Q) + A.U, 



with 



AM 



1.1 



ngK 2 ) 



(53) 



(54) 



where the symmetrization is necessary to achieve gauge invariance. Of course, the diagonal 
elements of AMjj correspond to the difference between the pole and running masses 
defined in eq. (|5T|). 
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6.2 Gauge Invariance 



As explained in section 2.2, the one-loop corrected vacuum expectation values are found 



by solving the one-loop corrected tadpole equations in eq. (|TT|) . Of course, it is desirable 
to work with gauge invariant vevs. In order to achieve the gauge invariance of the f Q 's, 
the one-loop tadpole T® R (Q) must be independent of the gauge parameter £. As it is 
shown in the appendix |C] the following set of tadpoles is gauge invariant: 




where S® denote neutral scalar bosons in the weak basis (see appendix A) 77's are the 
Fadeev-Popov ghosts. A similar set for the Z gauge boson exists. Nevertheless, the 
tadpole with a charged Goldstone boson in the loop introduces a gauge dependence that 
cannot be canceled. For this reason, the Goldstone boson loops are removed from the 
tadpoles T a (Q) and introduced into the self energies. This in turn allows us to achieve 
the gauge invariance for the two point functions, as well as for the vev's, as explained 
below. 

Among the loops contributing to the self energies, consider for example the VT-boson 
loop, which in the general gauge is 



W 



where the dots indicate terms proportional to 75 which are irrelevant for us, are 
charged fermions resulting from mixing between charginos and charged leptons, and 



X q—2 E (OffiOffi + OSgZOSS) |2S 1 (p 2 , ml m 2 w ) + B (p\ ml m 2 w ) 



-iB G { J) \mUml,)- r ^^- 



rn 



w 



Bi{p ,m W i m k) ~ B i(P ,t m wi m k] 



(55) 



(nS) W = ^ E {OlTkO^i + OfflOfii) m k [3B (p\ ml m 2 w ) + ^ (p 2 , ml ^m 2 w ) 

This graph introduces an explicit dependence on the gauge parameter £. The other self 
energy graph with £ dependence is the one that includes the charged Goldstone boson. 
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The charged Goldstone boson is one of the eight charged scalars S£ resulting from mixing 
between the two charged Higgs fields and the six charged sleptons. This contribution is 



S r 




IF- 



V\S+ 



+ 



where again, the dots indicate terms proportional to 75, and 

^ 8 5 
i07T r =l k=l 

( n 5) 5+ = -iiEEtW + W)" 1 *^^) (56) 

i07r r=l k=l 

with the couplings given in Appendix B. Nevertheless, gauge dependence is not canceled 
after combining eqs. (|55|) and (|56f ). In order to achieve it the inclusion of the Goldstone 
boson tadpole graphs, left over from the tadpole equations, is necessary: 



"1 " 



1.1 

) 

i . I . " i 



I CO 

^■0 1 °k p0 



where (Sj<) Tad = and 



K) Tad = -^J2{Ol^k + O^ k )^gil% s :A (^ w ) (57) 
6Z7T k=1 m s o 



J k 

,s°s+s- 



The A 0} Bq and B\ appearing above are Passarino-Veltman functions [31], g^Q+Q- being 
the neutral scalar coupling to a pair of charged Goldstone bosons, and O nns the neutral 
scalar couplings to a pair of neutral fermions (neutralino/neutrinos). Numerically, we 
have checked that, by adding the Goldstone tadpoles to the self energies our results do 
not change by varying the gauge parameter from £ = 1 to £ = 10 9 , thus establishing the 
gauge invariance of the calculation. Similarly we have also checked that the corresponding 
set of diagrams involving the neutral gauge boson tadpole (Z) plus neutral ghost tadpoles 
is gauge invariant and, similarly, the contribution to the self energies due to Z exchange 
plus neutral pseudoscalars and neutral Goldstones is also gauge invariant. 

Before we close this section we would like to add a short discussion on the basic 
structure of the loops which will be useful in the following. It is useful to do this in the 
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approximation where the ft p parameters are small, as discussed above. As seen from 
the expression for m e ff, at tree-level the effective neutrino mass matrix in this limit has 
the structure ~ AjAj, and at this level of sophistication neutrino angles are simple 
functions of ratios of Ai/ Aj. The one-loop corrections, however, in general destroy this 
simple picture. This can be seen as follows. The one-loop corrections have the general 
form, 

(Ey, Tkj) ~ YX°^°if> + 0i,<A,d)(Bi, mB ) (58) 
where the O stand symbolically for the various couplings. Now, since the expansion 



matrix £, defined in eq. ( p3[) can be written as ^ a ~ f a €i + 9aAi (see eq. (|35|) ) a product 
of two couplings involving neutrino-neutralino mixing has the general structure, 



O ha O hb ~ (f a ei + g a Ki) x (fa + g' a \ 3 ) x F{...) (59) 

where all the other dependence on the SUSY parameters has been hidden symbolically in 
F(...). The one-loop corrections therefore also carry a certain index structure, which can 
be written as 

m l - lo °p ^ ae . e . _|_ bfaAj + Ai€j) + cAiAj (60) 

where a, b, c are again complicated functions of SUSY parameters involving couplings, the 
Passarino-Veltman functions etc. Clearly, the terms proportional to c in eq. (0) above 
will lead only to a renormalization of the heaviest neutrino mass eigenstate. On the other 
hand the terms proportional to a in eq. (|60| ) are genuine loop corrections. Consider the 
simple case where all Aj = 0. Clearly in this case the tree- level neutrino mass is absent, 
but the one-loop effective neutrino mass has the same index structure as before, but now 
in terms of e^'s instead of Ajj's. In this idealized case angles are given as simple funtions 
of 6j ratios. For non-zero Aj the terms proportional to b in eq. eq. flSHD , however, destroy 
this simple picture. Any mismatch between e»/ej and Aj/Aj will lead, in general, to a 
very complex parameter dependence of the neutrino angles. 



7 Numerical Results on Neutrino Masses and Mix- 
ings 

Here we collect our numerical results on neutrino masses and mixings. As we have seen, 
a characteristic of the BRpV model is the appearence of vacuum expectation values for 
the sneutrino fields, v\ which imply a tree-level mass for one of the neutrinos given by eq. 



(p|). The one-loop-corrected neutrino mass matrix gives important contributions to the 
heaviest neutrino mass which we have determined through the renormalization procedure 
sketched above. 
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I A I 10 5 |A|/(VM2/i) 

Figure 1: Example of calculated AmJ ta as a function of (left) the alignment parameter A 
and (right), as function of |A|/ (^/M^/J,), all of these expressed in GeV. The figure shows 
that eq. (^Tj) can be used to fix the relative size of R-parity breaking parameters to obtain 
the correct AmL. 



First we note that in order to solve the atmospheric and solar neutrino problem one 
requires m 1 ~ loop -C m tree . If this is fullfilled it is essentially trivial within our model to 
solve the atmospheric neutrino problem. It is simply equivalent to choosing an adequate 
size of the alignment vector |A| 2 , as can be seen from eq. ([IT]) and is also demonstrated in 
Fig. ([I]). However there are regions of parameters where the one- loop contributions are 
comparable to the tree-level neutrino mass. This is discussed in quite some detail below, 
where we give an illustrative parameter study in order to isolate the main features of the 
dependence on the underlying parameters. First we get a rough idea of the magnitude of 
the neutrino masses including the one-loop corrections by displaying in Fig. @ the three 
lightest eigenvalues of the neutrino/neutralino mass matrix as a function of the parameter 
|e 2 |/|A|. Other parameters are fixed as follows: a) MSSM parameters: mo = fi = 500 
GeV, M 2 = 200 GeV, tan/? = 5, B = -A = m . b) RPV parameters: |A| = 0.16 GeV 2 , 
10 A e = A M = A T and e\ = e 2 = €3. In the left panel we give the predicted masses in the 
general case, while in the one on the right we apply the sign condition, 

(e M /e r ) x (A M /A T ) < (61) 

to be discussed in more detail below. 

One notices that the parameter |e 2 |/|A| determines the importance of the loop con- 
tribution relative to the tree-level-induced masses. For example, from the right panel one 
sees that, below |e 2 |/|A| <C 10 the heaviest neutrino mass m 3 is mainly a tree-level mass, 
while for |e 2 |/|A| ^ 10 the loop-induced masses are important relative to the tree-level 
one. Similar results are obtained for other choices of MSSM parameters. 
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Figure 2: Example of calculated neutrino masses in units of eV as a function of |e 2 |/|A|, 
for a particular though typical choice for the other parameters (see text), illustrating the 
relative importance of tree versus loop-induced neutrino masses. 

It is also interesting to analyse the dependence of the neutrino mass spectrum ob- 
tained in this model as a function of other supersymmetric parameters. In Fig. (|3|) we 
show the three lightest eigenvalues of the neutrino neutralino mass matrix as a function 
of tan/3, keeping the other parameters fixed as in Fig. (Q), fixing e 2 /|A| = 1. Again in 
the right figure we have applied the sign condition discussed in more details below. Loop 
contributions are very strongly correlated with tan (3. Similarly one can compute the three 
lightest eigenvalues of the neutrino/neutralino mass matrix as a function of m , as shown 
in Fig. (|j). Larger mo leads to smaller loop masses, as expected. From Fig. (0)-Fig. (^) 
one sees that, as expected, the pattern of neutrino masses obtained in the bilinear Jji v 
scenario, for almost all choices of parameters, is a hierarchical one. 

In the above we have not paid attention to whether or not the parameter values used 
in the evaluation of the neutrino mass spectrum are indeed solutions of the minimization 
tadpole conditions of the Higgs potential. We now move to a more careful study of the 
magnitude of the neutrino mass spectrum derived in the 1/ip scenario. 

In order to proceed further with the discussion of the solutions to the solar neutrino 
anomalies in this model we must distinguish two cases: 

1. unified universal boundary conditions on the soft SUSY breaking terms (SUGRA 
case, for short) 

2. non-universal boundary conditions on the soft SUSY breaking terms (MSSM case, 
for short) 

In what follows we refer to these two possibilities as SUGRA and MSSM cases, accordingly. 
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Figure 3: The neutrino mass spectrum versus tan /3, for parameters otherwise chosen as 
in Fig. (0). The importance of loops increases strongly with tan/3. 
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Figure 4: The neutrino mass spectrum versus mo, for parameters otherwise chosen as in 
Fig. (0). The importance of loops decreases with increasing m . 
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Figure 5: Am^ f versus |e|//x for \i < (left) and /i > (right). 



For the analysis of the neutrino masses these two scenarios are very similar so we 
focus on the case where the low-scale paramaters are derivable from a universal super- 
gravity scheme. In Fig. (|5|) we show the mass squared difference Am^ which is relevant 
for the analysis of neutrino oscillations and therefore relevant to the interpretation of solar 
data, as a function of the parameter \e\/ fi. In the left panel we display \i < while on the 
right panel \i > 0. Small values prefer Am^ in the range of the vacuum solution to the 
solar neutrino problem, while large values give masses in the range of the MSW solutions. 

Points shown in the following figures were obtained scanning the relevant parameters 
randomly over the region: M 2 and |/i| from to 500 GeV, m [0.2 TeV, 1.0 TeV], and 
bo [-3,3] and tan/5 [2.5,10], and for the Ift p parameters, |A M /A T | = 0.8 — 1.25, e M /e r = 
0.8 - 1.25, |A e /A r | = 0.05 - 0.1, e e /e T = 0.6 - 1.25 and |A| = 0.05 - 0.12 GeV 2 . They 
were subsequently tested for consistency with the minimization (tadpole) conditions of 
the Higgs potential and for phenomenological constraints from supersymmetric particle 
searches. 

One can also explicitly determine the attainable range of Am^ for which the corre- 
sponding A777I3 (see below) lies in the range required for the correct interpretation of the 
atmospheric neutrino data. The result obtained is displayed in Fig. (|(J) in which we show 
Am^ 2 as function of tanj3 for those points which solve the atmospheric neutrino problem. 



We now turn to the discussion of the three neutrino mixing angles and of how they 
must be identified in terms our our underlying parameters. Following the usual convention 
the relation 

M lL = W M (62) 
v a = U ak v k (63) 
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Figure 6: Am^ versus tan/3 for points which solve the atmospheric neutrino problem. 



connecting mass-eigenstate and weak-eigenstate neutrinos are recovered in our notation 

as 

= MgkMa (64) 

where the mixing coefficients M are determined numerically by diagonalizing the neutral 
fermion mass matrix. Note that, without loss of generality, in the bilinear model one can 
always choose as basis the one in which the charged lepton mass matrix is already diagonal. 
The neutrino mixing angles relevant in the interpretation of solar and atmospheric data 
are identified as (if U e $ C l,as indicated by the atmospheric data and the reactor neutrino 
constraints). 

sin 2 (20 13 )=4[/ 2 (l-[/ 2 ) (65) 



sin 2 (20 12 ) = AU^U, 



e2 



(66) 



Note that the maximality of the atmospheric angle is achieved for A M = A r (see Fig. ([/])) 
and A e is smaller than the other two, as required by the Chooz data (see below). In fact 
we have found f32fl , that if e 2 /A <C 10 then the approximate formula holds 



U 



o3 



A Q /|A| 



(67) 



In Fig. (|D we show the expected magnutide of U^ 3 versus the relevant ratio of Ifi p 
parameters. In order to comply with the reactor data from the Chooz experiment one 
should have ?7 2 3 below 0.05. This implies a bound on A e which can be read off from the 
figure. 

The discussion on the solar mixing angle is more involved. First note that it has no 
meaning before adding the one-loop corrections to the neutrino mass, since in that limit 
the two low-lying neutrinos would be degenerate in mass. 
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Figure 7: Atmospheric angle versus A^/yA^ + A^. Maximality is obtained for A M ~ A T 
if A e is smaller than the other two (see Fig. (§)). 
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Fi gure 8: t/g 3 versus A e /yA^ + A^. To obey the experimental bound < 0.05 A e must 
be smaller than A M , A T . 
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Figure 9: sin 2 (28 so i) versus A e /JA 2 + A 2 for the SUGRA case, for a discussion see text. 



In order to proceed further with the discussion of the solutions to the solar neutrino 
problem in this model we must analyse carefully the implications of eq. (|2"B"|). Here it is 
important to distinguish between case 1 (SUGRA) and case 2 (MSSM) discussed above. 



In the SUGRA case by taking the ratio of the first two equations in eq. 

e e Am 2 — tan /3[iAB e A e 
e„ Am 2 - tan/3/iAB^ A, 



(68) 

we conclude that, since A e < A M and, since the relevant ratio of SUSY Soft-breaking 
terms is close to one, it follows that sin 2 (2# ) is small. The predictions for the solar angle 
as a function of the lji p breaking parameters is indicated in Fig. @. 

More precisely, the interpretation of the solar data || in terms of the small angle 
MSW solution indicates that 

sin 2 (2fl ) < 1(T 3 - KT 2 (69) 

and this in turn selects the required ratio of A e to A M and A r . Therefore in this case the 
large angle solutions, including the vacuum or just-so solutions do not fit in the scheme. 

We now move to the general MSSM case. In this case the ratio of SUSY soft- 
breaking terms appearing in eq. ([BSD is in general arbitrary and thus the ratios of Aj/Aj 
is no longer tied up to the ratios of ei/e/s. This opens up the possibility for large angle 
solutions to the solar neutrino problem. At first sight it would seem that all predictivity 



of the solar angle is lost in this case, as seen in left panel of Fig. ([10) . 



The ability of our model to determine the solar neutrino angle may be understood 
in terms of eq. (pi). For example in the SUGRA case we see from eq. fl68|) that the e 
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Figure 10: sin 2 (28 so i) versus e e /ye^-Tef . The left panel corresponds to the case without 
the sign condition and the the right panel assumes the sign condition. 

and A ratios are fixed within a narrow range, leading to the small mixing angle prediction 
for the solution to the solar neutrino problems. There is however another way to obtain 
predictivity for the general MSSM case, namely by applying eq. (|6T|). 

The possibility of our model predicting the solar angle even in the general MSSM 
case by assuming eq. (|6lD can be understood as follows. Consider first the simplified limit 
A e = 0. In this case v\ = v e at tree-level and there is no mixing at all between the electron 
neutrino and the other two states, but a finite mixing exists at one-loop, due to the terms 
proportional to e e . In this case the sign condition, defined in eq. (|6l| ) introduces two more 
zeros into the matrix proportional to b in eq. fl60"|) above, if |e M | = |e r | and \A^\ = |A T |. 
This fact simplifies the calculation of the solar angle very much, since one of the neutrino 
eigenvectors (the one for v e ) has no dependence on the Aj ratios but only on the e, ratios. 
For a non-zero A e (and small departures from equality of e M ,e T or |A M |, |A r |) this feature 
is destroyed and a A e dependence reintroduced in the solar angle. However, as long as the 
one-loop contributions are smaller than the tree-level one and as long as A e <C A^ T , the 
"cross-talk" between the A e and e e pieces is sufficiently small, such that some predictivity 
of the solar angle is retained, as illustrated in figure Fig. fliCf ) (right panel). 

The discussion on mixing angles may be summarized as follows. In the case that 
one-loop corrections are not larger than the tree-level contributions, the approximate 
formula 

U a3 « A Q /|A| (70) 

holds. This allows one to fix the atmospheric angle and at the same time obey the CHOOZ 
constraint. For the solar angle, however, the results depend on whether one wants to work 
in a SUGRA motivated scenario or not. For the SUGRA scenario we have found that 
our model allows only the small mixing angle MSW solution (SMA), while for the general 
case also LMA and vacuum oscilation solutions are possible. 
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8 Conclusions 



We have shown that the simplest unified extension of the Minimal Supersymmetric Stan- 
dard Model with bi-linear R-Parity violation typically predicts a hierarchical neutrino 
mass spectrum, offering a natural theory for the solar and atmospheric neutrino anoma- 
lies. In this model only one neutrino acquires mass due to mixing with neutralinos, while 
the other two get mass only as a result of radiative corrections. We have performed a 
full one-loop calculation of the effective neutrino mass matrix in the bi-linear lfi p MSSM, 
taking special care to achieve a manifestly gauge invariant calculation and performing 
the renormalization of the heaviest neutrino, needed in order to get reliable results. The 
atmospheric mass scale and maximal mixing angle arise from tree-level physics, while the 
solar neutrino scale and oscillations follow from calculable one-loop corrections. 

Under the assumption of universal boundary conditions for the soft-supersymmetry 
breaking terms at the unification scale we find that the atmospheric scale is calculable by 
the renormalization group evolution due to the non-zero bottom quark Yukawa coupling. 
In this case one predicts the small mixing angle (SMA) MSW solution to be the only 
viable solution to the solar neutrino problem. 

In contrast, for the general MSSM model, where the above assumptions are relaxed, 
one can implement a bi-maximal [[33] neutrino mixing scheme, in which the solar neutrino 
problem is accounted for through large mixing angle solutions, either MSW or just-so. 
A great advantage of our approach is that the parameters required in order to solve the 
neutrino anomalies can be independently tested at high energy accelerators, as originally 



proposed in ||14|| . In fact, as shown in ref. Jj2[ |3J] the bilinear ^? p model predicts the lightest 



supersymmetic particle (LSP) decay to be observable at high-energy colliders, since the 
expected decay path can easily be shorter the typical detector sizes. This happens despite 
the smallness of neutrino masses indicated by the SuperKamiokande data. This provides 
a way to test this solution of the atmospheric and solar neutrino anomalies and potentially 
discriminate between the large and small mixing solutions to the solar neutrino problem. 
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A Mass Matrices 



A.l Scalar Mass Matrices 
A. 1.1 Charged Scalars 

The mass matrix of the charged scalar sector follows from the quadratic terms in the 
scalar potential 



V guadratic = S'-Ml±S' + (71) 

where the unrotated charged scalars are S' + = (H£, H+, ej, fi'^f^, J^rTr)- For conve- 
nience we will divide this (8 x 8) matrix into blocks in the following way: 



Ml 



M 2 HH tfgf 



M 



Hi 



+ t,rn 2 w 



M\ M 2 J 
M 2 B M 2 C 



(72) 



where the charged Higgs block is 



M 2 HH 



Bfj, + \g 2 v d v v 



+^El 1 e i f d + lEl j =iV i (h E h E ) 



■'j 



Bfi + \g 2 v d v u 



(73) 



B^+y'ivj+EUv 2 ) 

E3 p , Vj i t u 
i=l J -'i t i „ ~r „ 

This matrix reduces to the usual charged Higgs mass matrix in the MSSM when we set 
Vi — €i — and we call m\ 2 = B^. The slepton block is given by 



m \l M 2 lr - 
Mrl M 2 rr1 



(74) 



where 
(M 



LL 



H (h E h T E ).. + \g 2 (- J2 vl-v 2 d + vlj Sn + \g 

+fi—ei5ij - a ( V] —e k I 5ij + aej + M 2 Lji 
Vi \ k =i^ J 

; i "i 



" ViVj Bieida H 6 



'■j 



M'ui = (v d A* E - fiv u h* E ) 
M 2 RL =(M 2 LR f 



(75) 

(76) 
(77) 



(M^).. = W 2 {-p4-v 2 d + v^5 ij + lv 2 d (hlh E 
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\fc=l / \s=l 



+ M 2 Rji 



(78) 



We recover the usual stau mass matrix again by replacing v j = 6j = (note that we need 
to replace the expression of the tadpole t{ in eq. @ before taking the limit). The mixing 
between the charged Higgs sector and the slepton sector is given by the following 6x2 
block (repeated indices are not summed unless an explicit sum appears): 



-fj£i - ±V d 



( h *E h E) ik v k + \g 2 v d Vi 



k=\ 



Vk 



-B^ + \g 2 v u Vi 



(-kVd) 



k=l 



(79) 



k=l k=l 

and as expected, this mixing vanishes in the limit = = 0. The charged scalar mass 
matrix in eq. (|72|), after setting t u = t d = U = 0, has determinant equal to zero for 
£ = 0, since one of the eigenvectors corresponds to the charged Goldstone boson with zero 
eigenvalue. 

For our one loop calculations one has to had the gauge fixing. The part of the mass 
matrix in Eq. (|72|) that comes from the gauge fixing reads for the (2 x 2) A block 





r vj 


-v u v d l 


M\ = 


V 2 


V 2 




-v u v d 


v 2 

u 




V 2 


V 2 



10) 



for the (6 x 2) B and the (6 x 6) C blocks 

ViV d -ViV u 

r2 



M 



B 







v 




M 



c 



M 



D 



(81) 



where the (3 x 3) D block is 



M 



D 



r i 


ViV 2 


VlV 3 -| 


V 2 


V 2 


V 2 


V 2 Vi 


vj 


V2V3 


V 2 


V 2 


V 2 




V2V3 


vj _ 


V 2 


V 2 


V 2 



(82) 



The charged scalar mass matrices are diagonalized by the following rotation matrices, 

,2 _2 N dS± 



with the eigenvalues diag(mg l , . . . , m| g ) 



i2 s Ml, 



(« S ) 



A. 1.2 CP-Even Neutral Scalars 



The quadratic scalar potential includes 



quadratic 



,0 ~R 



V; 



}M 



S° 



On 



+ 



54) 
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where the neutral CP-even scalar sector mass matrix in eq. fl84|) is given by 



M 



s° 



M 



ss 



M 



SVR 



where 



M 2 SS 



19 9 x — ^ 

B^— + jg z v d + 11 2^ e fc— + — 



-Bfi - \g 2 z v d v u 



M 



-fiet + \g 2 z v d Vi 
Bi€i - \g 2 z v u Vi 



(86) 



*7) 



and 
v 



2 _ 



'J 



-- fi€i Biti €i ^ e k | E _ ( M «fc + M , 



where we have defined g\ = g 2 + g' 2 - In the upper-left 2x2 block, in the limit t> j = = 0, 
the reader can recognize the MSSM mass matrix corresponding to the CP-even neutral 
Higgs sector. To define the rotation matrices let us define the unrotated fields by 



S 



( „0 ~R ~R ~R\ 



Then the mass eigenstates are S® given by 



CO _ t>S° q/O 



with the eigenvalues diag(mg l , . . . , m| 5 ) = R s ° M\ (-R S °) • 



A. 1.3 CP-Odd Neutral Scalars 



(89) 



(90) 



The quadratic scalar potential includes 



V r , 



quadratic 



+ 



where the CP-odd neutral scalar mass matrix is 



M 2 p0 



M 
M 



pp 

Puj 



M 
M 



Puj 

2 _ 

VlVl 



Ml 
M 2 F 



M 
M 



2T 
2 



(91) 



(92) 



where 
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and 



M 



pp 



Bp h p > e k 1 

v d ^[ v d v d 

Bp 



Bp 



Bp 2^ B k e k 1 

Vu k=l ^ u ^ u 



-pti 
— I>,< ; 



(93) 



(94) 



Ate* -Di^i £j 2^ e fc — 



2 2^ 

fe=l y * 



^ (Mf tt + M 

11; V 



2 

Lki 



+e . e . + i ( M |.. + M ; 



2 



(95) 



Finally the part of the mass matrix in Eq. (|9~2"D that comes from the gauge fixing 
reads for the (2 x 2) E block 



M 



v 



-v u v d 



for the (3 x 2) F block 
and for the (3 x 3) G block 



Mi 



V 2 V 2 
ViV d -ViV u 



(96) 



(97) 



M 



G 





VlV 2 




V 2 


V 2 


V 2 


V 2 Vi 


vj 


V2V3 


V 2 


V 2 


V 2 




V3V2 


vj _ 


V 2 


V 2 


V 2 



(9* 



The charged pseudo-scalar mass matrices are diagonalized by the following rotation ma- 
trices, 

(99) 

,o\T 



with the eigenvalues diag(m| l , . . . , m| 5 ) = R p ° M 2 p0 (R P °) . where the unrotated fields 

p'° = {vl vl 2?, (100) 



arc 



A. 1.4 Squark Mass Matrices 

In the unrotated basis u[ = (uLi,u* Ri ) and d\ = (d,Li,d* Ri ) we get 

Viatic = \u' ] M~ 2 v! + \& M~ d! 



(101) 
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where 

(MX, Ml \ 
2= « Li M (102) 

with q = (u, d). The blocks are different for up and down type squarks. We have 

M Ll = \ v l h *uh T u + Mq + | (Am 2 w - m 2 z ) cos 2(3 
M Zrr = \ < huhh + M 2 + 2 3 (m 2 z - m 2 w ) cos 2(3 

MIhl = MIlr (103) 



and 



M Sll = ~vlh* D hl + M 2 Q -l(2m 2 w + m 2 z )cos2(3 
M Im = \ v > T D h D + M 2 D - \{m 2 z - m 2 w ) cos 2(3 
M~ = ^=A* D -^h 

dLR ,/9 u 1 

t 



V2 D ^V2' 



D 



dRL dLR 

We define the mass eigenstates 
which implies 



Ml BT = M1J (104) 

q = fflq i (105) 



& = R q 3i q (106) 
The rotation matrices are obtained from 

W ] (Mf 19 ) 2 W = M~ 2 (107) 

In our case the matrices in Eq. (|102| ) are real and therefore the rotation matrices R q are 
orthogonal matrices. 



A. 2 Chargino Mass Matrix 

The charginos mix with the charged leptons forming a set of five charged fermions Ff 1 , i = 
1, . . . , 5 in two component spinor notation. In a basis where i( +T = e^, /i^>, Tr) 

and i(~ T = (— i\~ , H~[ , e~[, /j,~£, t£), the charged fermion mass terms in the Lagrangian are 

where the chargino/lepton mass matrix is given by 
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Af, 



c 



' M 


-h=QV, t 













f 


1 / 7 \ 


1 / 7 \ 

-73(^)22^2 


1 / 7 \ 

-75(^)33 ^ 




-ei 










T2 gV2 







73(^)22 v d 







-63 








73(^)33^ 



(109) 



and M is the SU{2) gaugino soft mass. We note that chargino sector decouples from the 
lepton sector in the limit 6j = — 0. As in the MSSM, the chargino mass matrix is 
diagonalized by two rotation matrices U and V defined by 



Then 



U f M c V 1 = M C D 



(110) 



[111) 



where Men is the diagonal charged fermion mass matrix. To determine U and V we note 
that 

Mq D = VM^McV 1 = WMcM^W) 1 (112) 

implying that V diagonalizes MqMc and U* diagonalizes McMq. For future reference 
we note that 

^-=U* kj F k ~ ; ^ = V kj F+ (113) 

In the previous expressions the are two component spinors. We construct the four 
component Dirac spinors out of the two component spinors with the conventions"'", 



Xi 



(114) 



B The Couplings 

B.l The Neutralino Couplings 



Using four component spinor notation the relevant part of the Lagrangian can be written 

as 



c = xi Y (o%[P L + o™p r ) x °j w~ + Y (0% p l + o%[P R ) x] w+ 

tHere we depart from the conventions of ref. |2f| because we want the e~ , fx~ and t~ to be the 
particles and not the anti-particles. 
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+Xi {0 C Lf jk PL + O c ^ jk P R ) x°j s k + Xi (pij ik PL + Olf ik P R ) Xj s£ 
(O^Pl + O^Pr) X° Z° + £J (0^ k P L + O^Pfl) x? Hi 

(02^ fc PL + OS5 fc P fl ) x? ft + xF (o^Pl + og fc P*) <& £ 
where g can be either d or u. The various couplings are: 

B.l.l Chargino— Neutralino— W 

/ 3 



o% = grans 



k=i 



nncw / /OcnwA * . /Oncw / /OcnwA 

U Lij - [ U Lji J ' ^ffij - ^iiji J 



B.1.2 Neutralino-Neutralino-Z 



B.1.3 Chargino— Neutralino— Charged Scalar 

/Qcns _ 



^ (h E11 N* j5 V* i3 + h E22 N* j6 V* i4 + /isasJV^Vis) 



— J2f 3 iiV* ( ;iV*,;{ — i?f 4 I1E22N* j 3 V* i4 — R k5 h,E3zN* j 3 V* i5 
- Rlt g' s/2N* 3l V\ 3 - R% g' s/2N* n V* lA - R s k t 9 'y/2N* n V* i5 
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/ncns _ 

u Rijk — <h 



+ R s k t (^ 2 N j2 U i3 + ^N iX U a - gN j5 U a ) 
+ R s k t {^N j2 U i4 + ^N^U* - gN j6 U a ) 
+ R s k t (^ 2 N j2 U i5 + ^N jX U a - gN j7 U«) 
+ R s k th E11 (N j5 U i2 - N j3 U i3 ) + R s k7 h E22 (N j6 U i2 - N j3 U l4 ) 
+ R S k th E 33 (N j7 U t2 - N j3 U i5 ) 



r~\ncs / ^ncs \ * . s^icns / ^ncs \ " 

U Lijk — [yRjik J i U Rijk — y^Ljik J 



B . 1 .4 Neutralino Neutralino— Scalar 



/nnnh _ _ 
U Lijk — 'Ij 9 



i2g (-gN* i2 N* j3 + g'N* a N* j3 -gN* j2 N* l3 + g'N* n N\ 3 ) 



+ J*£ (+ 5 JV* i2 JV* j4 - </ JVixJV*^ + 5 N* j2 N* i4 - g' fiT^N**) 
+ R s k3 (-g N* l2 N* j5 + g' N* a N* j5 - g N* j2 N* i5 + g' N* j± N* i5 ) 



+ R s k4 (-g N* t2 N* j6 + g' N* n N* j6 - g N*, 2 N% 6 + g' N* n N* i6 ) 



+ Rf 5 (- 9 N% 2 N* 37 + g'N* a N* j7 - g N* j2 N* i7 + g'N* n N* i7 ) 



nnnh ( /^)nnh "\ 

•^ffijfc — \yLjik ) 



B.1.5 Neutralino-Neutralino-Pseudo Scalar 



rinna 
U Lijk 



R p kl (-gN* i2 N* j3 + g'N*^*^ -gN* j2 N* l3 + g'N*^*^ 



+ iC (+gN* l2 N* jA - g'N* a N* j4 + gN* j2 N* l4 - g , N* jl N*^) 



+ iC (-gN* l2 N*, 5 + g'N* a N* j5 - g N* j2 N\ b + g'N* n N\ 5 ) 



+ iC (-gN* t2 N* j6 + g'N* a N* j6 -gN* j2 N* l6 + g'N*^*^ 



+ iC (- g N* i2 N* j7 + g'N* a N* j7 - 9 N* j2 N\ 7 + g'N* 3l N* a ) 



/Qnna / nnna \ " 

^Rr/fc — [yijik ) 
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The factors rji are the signs one has to include if we consider N, U and V as real matrices 
and the mass of the fermion i is negative. 



B.1.6 Neutralino-Up Quark-Up Squark 

°lTjk = 3(^5) tan() wN* ' jlRk,m+3^R itTn ~ (K)mlRk,m R Rij N *j4 (123) 

(K) ml RT +3il R* u Lm .N j4 (124) 

{O^Y (125) 
B.1.7 Neutralino— Down Quark— Down Squark 

OlT jk = tan^V< m+3 ^ m - (h d ) ml Rl m R Ri ^* j3 (126) 



02% = -(fy( N * + 3 tan^O^^L ti ~ 

and, 

^ujk — [ynjik ) i w Rijk — 



Ot 3 k = (^)(^ --t a n6 w N n )Rt tm R* d Lmi - (h:) ml 4] +3 R* d Lm .N j3 (127) 

and, 



B.2 The Neutral Scalar Couplings 

To evaluate the tadpoles we need the couplings of the neutral scalars with all the fields in 
the model. These couplings are easier to write in the unrotated basis. The couplings for 
the mass eigenstates can always be obtained by appropriate multiplication by the rotation 
matrices. As an example, and to fix the notation (repeated indices are understood to be 
summed unless otherwise stated), the couplings of three neutral scalars in the two basis 
will be related by 

n S°S°S° nS° pS° TfS° n S'°S'°S'° /-, 9Q ^ 

9ijk ~ K ip K jq K kr 9pqr { iZ ^> 

Sometimes we will also use partially rotated couplings, for instance 

n S' S°S _ rfS { > pS° n S ,0 S ,0 S'° /-. on\ 

9ijk ~ 11 jq "'kr 9iqr \ LOV J 

in an obvious notation. These couplings are defined as follows 

s°s°s° d^jC 
9l3k = dSfdSfdSf 

s ,0 s'°s' d 3 £ 
9ijk - dS'°dS'°dS' k ° { } 
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B.2.1 Neutral Scalar-Neutral Scalar-Neutral Scalar 



9ijk S '° s ' = -\ {p 2 + a' 2 ) u m (5 m % Sjk + $mj s ik + s mk 4,) 



where we have defined 

u m = (y d ,v v ,v 1 ,v 2 ,v 3 ) 
For future reference we also define 



Sij = diag(+, -, +, +, +) 



V m = (vi,V 2 ,V 3 ) 



while 8^ without the hat is the usual Kronecker delta. 



B.2.2 Scalar-Pseudo Scalar-Pseudo Scalar 



(132) 



(133) 



(134) 



9ijk = ~4 [9 +9 ) u m 5 mi 5 jk 



(135) 



B.2.3 Scalar-Charged Scalar-Charged Scalar 



We define 



/ n S'°S'+S' 
/ 9iHH 



n S'°S'+S' 
9ij k 



n S'°S'+S'- 
9iHL 



n S'°S'+S'- v 
9iHR \ 



n S ,0 S'+S'- 
HLR 



f n s ,0 s'+s'-\" i n S'°S'+S'- 

\9iHL ) 9iLL 9, 

i ^5'°5'+5'-V „S'°S"+S"- ; 

V \9iHR ) \9iLR ) 9iRR ' 



where 



f„s'°s'+s'-\ 

\9iHH ) 



\9 2 



(136) 



- v u (8u Sji S k2 + S i2 Sji 5 kl + 8u S j2 5 kl + 5 i2 S j2 5 k2 ) 

-Vd ($n Sji 5 k i + Sn 5 j2 5 k2 + S i2 5 jX 5 k2 + 8 i2 5 j2 5 kl ) 
+ v m 5j_ 2)Tn (Sji S k i — Sj 2 S k2 ) 

l /2 ? c 

4 9 r ^ j m < J%m < Jj k 



-\v m SjiS k i (h E h E + h* E h 



m,i—2 



(137) 



/ 's m s'+s'-\ 

\9iHL ) . 



is 2 



Sij S mk ~\- 2 V m (h E h T E ) mk SuSji 



+\v d Sji (h* E h E ^._ 



i-2,k 



(138) 



f„S'°S'+S'-\ 
\9iHR ) 



, = 73 (^) mfc ^2 + <^i) + 73 (^E)i- 2 , fe ^i 

+73^(%- 2)fe ^2 (139) 
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{9iLL +S ' ) jk = \ ~ 9' 2 ) U m 5 im 5 jk ~\g 2 V m (<Ji_ 2 J Smk + 6i- 2 ,k 

- {h* E h T E). k v d 8 a 
{£T S ") jk = -jsS a (A%) jk + ^„S a {h E ) 3k 
(9iRR +S ' ) jk = \g 2 u m l im b jk -v d b iX (h T E h* E )), k 

( h *E)i-2,k ( h E) mj + ( h E)i-2,j ( h * E )mk 



2 v m 



B.2.4 Scalar-Up Squarks-Up Squarks 



With the definition 



we get 



where 



S'°u'u'* 
9iLL 

n S'°u'u'* 
9iLR 

9iRL 

„S'°u'u'* 
9iRR 



£ — 9ijk U U u 'j u 'k* + 



n S'°Z'u'* 

9ijk 



n S'°u'u" n S'°u'u" 

9iLL 9iLR 

n S'°u'u'* n S'°Z'Z'* 

, UiRL 9iRR 



u m 5 im (-| g 2 + jz g' 2 ) I-v u (huh\j) 5 i2 
~ 73 ^ u ~*~ 73 V ^ u ^ a ~ 73 ^ u e ™ ^-2,™ 



n H'°u'u" 
9iLR 



~\u m lira g' 2 1 ~ V u (hjjh*^ 5 i2 

where X is the unit 3x3 matrix. 



B.2.5 Scalar-Down Squarks-Down Squarks 



With the definition 



we get 



where 



„S'°d'd'* 

giLL 

„S'°d'd'* 
UiLR 

n S'°d'd'* 
tliRL 

„S'°d'd'* 
fJiRR 



£■ — 9ijk d d $1° d'j d k + 



S'°d'd'* S'°d'd' 

S'°d'd'* I 9iLL 9iLR 

I n S'°d'd'* n S'°d'd" 

> UiRL 9iRR 



U m km (\g 2 + J2 9 12 ) 1-V d {h D hty S H 

A D + -^fih D 5 i2 



n S'°d'd>' 
9iLR 



\ u m 5 im g' 2 l - v d {h T D h*j^ 5 a 
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B.2.6 Scalar-W+-W- 
With the definition 

C = gf w+w ~ S? W + W~ + ■ ■ ■ (149) 

we get 

g s ,0 w+w = g ^ + Vu + Vm s._ 2 m } (150) 

where 

v = \jv\ + v\ + v\ + v\ + v\ (151) 

B.2.7 Scalar-Z°-Z° 
With the definition 

C=\gf W S'°Z°Z° + --- (152) 

we get 

ft = 7 (Ud<>il +W«O i2 + W m (3i_2, m ) (153) 

cos 9 W v 
B.2.8 Scalar-Quark-Quark 

With the definition 

C = g*T u Sf u 3 u k + 4? Sf d 3 4 + • • • (154) 

we get 

9i* U = "73 So (155) 

and 

4? = "75 (^)i* ^ ( 156 ) 

B.2.9 Scalar— Chargino— Chargino and Scalar-Neutralino— Neutralino 

With the definition 

£ = XT (OS! + Ogi) xj S<° + ^ XF (OZg£ + 0$ k ) x? S? (157) 

we have 

g {y*iiU* j2 4i + V* i2 U*ji 5 k2 + V",iZ7' /;i <5 fe3 



Lijk 



V2 



+ V*nU* j4: 5 k A + V*iit/*j 5 <5 fc5y 

+ (/l^llf *j3V*i3 + h E 7!lU* 34V* i\ + h E33 U* j5 V* i5 ) 5 kl 
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(h E11 U* j2 V*i3 5 k3 + h E22 U* j2 V* iA 5 k 4 + h.E;uU* j2 V* i5 8 kb ) 



^cch' 
w Rijk 



(p C Ljik) 



(158) 



and 



Oltk = V, g (-9 X* i2 N* j3 + g' N* tl N*, 3 - g N*, 2 N* t3 + g' N* n N* l3 ) 

(4l - S k 2 + 43 + ^4 + Sk5) 

{oi%y 



,onnh' 
W Rijk 



(159) 



C Tadpoles 

C.l Gauge Boson and Ghost Tadpoles 

We will consider the gauge boson and ghost tadpoles in an arbitrary gauge to show 
that the dependence on £ cancels out. We will do it for any model. 

C.l.l General Z° Boson Tadpole 

We write down the tadpole contribution from the Z° for a general theory with the coupling 
gnzz to the higgs boson. 




IgHZZ 

H 



iT z = \ig H zz J J^ G AP) 



(160) 



where 



<V = H) 



(i-0- 



\f-M\ ^ -'{jt-MlW-tMD (161) 
and the factor | is a symmetry factor. Now we do some transformations in the second 
term of the Z° propagator 

.2 1 



c(0 



;i-o 



1 



and therefore we can write 



-0 



3 1 

p 2 — Mf p 2 — £Mf 



(162) 



(163) 
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Then 



16tt 2 



3A)(Mf) + e M£M 2 Z ) 



16tt 2 2 

where we have used the definition 



9hzz 



3A (M 2 )+^A ^M 2 ) 



d d p 



16tt 2 



) d p 2 — m 2 



(164) 



(165) 



As A (£m 2 ) grows for large £ as £m 2 we conclude that grows like £ 2 . This dependence 
has to cancel against other diagrams. It is easy to realize that the Goldstone of the Z° 
will not do it because, although its mass depend on £, its contribution to the tadpole will 
only grow like £ because its coupling to H does not depend on £. But the ghost coupling 
to H does depend £ as we will see. 



C.1.2 General Z° Ghost Tadpole 



Let us then calculate the tadpole of the ghost of the Z°. We have 



CJ 



••t* 



igHC z C z 



iT Cz = (-l)i9H Cz c z J 



d d p i 



(166) 



where the factor (—1) is because of the anti-commutative properties of the ghosts. Using 
the definition of A we get 



^ = -^- 2 9Hc-c z A (iMl) (167) 



Adding the two contributions together we obtain 

% 



-g HZZ A (M z ) + [\g H zzi + 9h Cz c z ) A (^M 2 z ) 



16tt 2 

We see that for the £ dependence to cancel one must have 



(168) 



19hzz^ + 9hc-c z =0 (169) 

As we will show below this is true for the SM, MSSM and also for the Bilinear R-Parity 
Model. Then the contribution from the Z° and neutral ghost tadpoles is, for any model, 
gauge independent and given by 

iT z + iT Cz = -±-jl gHZZ A (Mj ! ) (170) 
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C.1.3 The W ± Boson and cF Ghost Tadpoles 



The calculation for the W boson and charged ghosts is very similar. The main differences 
are that the W tadpole does not have a factor | and that there are two ghosts for the 
W ± . Therefore we have 



iT w + iT r + + iT. 



w c w IQti 



3gnww A (M W ) 



(171) 



+ (gnww £ + 9 Hc +^ w + 9 Hc ^- w ) MiK 

We see that the £ dependence will cancel out if 

9hww £ + 9 Hc+w ~ w + 9 Hc - w ~ w = (172) 

We will show below that this is true in general. Then the contribution from the W ± and 
charged ghost tadpoles is, for any model, gauge independent and given by 

iT w + iT c+w + zT c - = 3 g HWW A Q {M 2 W ) (173) 
C.1.4 The Standard Model 

Now lets us see us the cancellation occurs in the Standard Model (SM). The relevant 
couplings for the Z° are 

9uzz = — ^— M z 
cos t>w 

9hc- Cz = ~-\-tM z (174) 
2 cos Ow 

and we immediatly see that Eq. ( |169| ) is verified. For the W ± we have 



9hww = gM] 



w 



9hc+ w ~ w = -\tM w 



9h c ,~ w = ~\ZM W (175) 



satisfying Eq. dl72|) . 



C.1.5 Bilinear R-Parity Model 

In the bilinear R-parity model the relevant couplings are 



"0^^~ y Vddil + Vu ° i2 + V m0i-2,m) 
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2 v 



and 



9 



9 



m z 



(176) 



(177) 



2 cos 9w ' v 

Then using Eqs. ( |1 50| . |1 58| . |1 76| . |1 77[ ) in Eqs.( |172| , |169| ) we see that the same cancellation 
occurs. 



C.2 General Tadpole Expressions 

After showing the gauge invariance of the gauge boson tadpoles together with their ghosts 
we give now the general tadpole in a compact form. We will write them for the unrotated 
neutral Higgs H'° because that is what is needed for substitution into Eq.(§). The general 
form can be written as (X = W ± , Z°, S ± , H°, A , u, d, u, d), 



where 



» ibit 



p w = 3gr w ^ W -A (M^) 

pf= l 9 r z ° zo MMi) 

8 

— - 2^ 9ikk A o{ m k) 

k=l 

P s ° - V- n s ' 0s0s0 AJm 2 ) 
k=i z 

5 i 

P L n S'°P°P a (2\ 

n - - 2^ T,9ikk A o{ m k) 
k=l z 

k=l 
k=l 

5 /Occh' I /~)cch' 

Pf = - E Lkki 2 Rkki 

k=l 



k=i / z 
^ = -E(-3)^"4m fc AoK) 



fc=i 

^ = - E (-3) 4m fc A (m 2 fc ) (179) 
fe=i 

where ' means that we sum over all fields except for the goldstone boson. As explained 
in section [O] the contribution of the goldstones is added to the self-energies to achieve 
gauge invariance. 
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D One Loop Self— Energies 

In this section we write down the contribution of the several self energy diagrams in the 
f = 1 gauge. 



D.l The W and Z Loops 

The contribution of the W and Z loops to the functions H v and U v can be written in 
the form (X = W,Z), 



y v 



n 



Y J F^ k B l {p\mlm 2 x ) 



16tt 2 V 13 

k 



16tt 2 



m 



x) 



;iso) 



with 



771 W 9 / nncw/ncnw , /nncwncnw\ 

r ijk — ^ \y Ljk w Lki ^~ ^ Rjk^ Rki ) 

G% k = -4 (0^ k O c R Z + OfftO™) (181) 

and 



^ijfc — \ U Ljk W Lki "T U Rjk U Rki ) 

Gf jk = -^(o^kO^ + O^OlZ) (182) 



D.2 The Scalar Loops 

All the scalar contributions can be written in the form (X = S ± , S°, P°, u, d), 



4 = -^EE4%>^ r 2 ) 

ng = -iEE G i m ^o(p 2 ,^,m?) (183) 

i07r r fc 

with 

p5 ± / ^ncs ^cns , ^ncs ^cns \ 

r ijkr — \y Rjkr u Lkir ~t~ ^Ljkr^Rkir ) 

— fD ncs D cns J- D ncs D cns \ ho/\ 

^ijfcr — \y Ljkr U Lkir U Rjkr U Rkir ) K 10 ^) 
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nnh /onnh 

Rkir 



1 / /^)nnh rinnh , /^jnnh /^) 

2 y^Ljkr^Lkir u Rjkr^ 

1 / /^jnnh /nnnh , /^\nnh 

2 \ W Ljkr U Rkir U Rjkr W 

1 / nima nnna , rimia /--) 

2 \y Rjkr U Lkir U Ljkr U 

1 / rinna nnna , rimia 

2 \yLjkr u Lkir u Rjkr w 



nnh 



nna 
ijfcir 



nna 
ijfcir 



uns 
ijfcir 



/ /onus ^uns , ^nus 
\yRjkr u Lkir w Ljkr u 

/^jnus z^juns . /^nus /'-juris 
W Ljkr^Lkir ~t~ ^ Rjkr W Rkir 



(^-jnds ^-jdns , A-jnds .odns 

^Rjkr^Lkir U Ljkr U Rkir 

( /^jnds /^jdns , /^nds /^dns 

\ U Ljkr U Lkir ^Rjkr^Rkir 
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